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ABSTRACT 

In this paper, the first of a series, we study the stellar dynamical and evolutionary processes 
leading to the formation of compact binaries containing white dwarfs in dense globular clus- 
ters. We examine the processes leading to the creation of X-ray binaries such as cataclysmic 
variables and AM CVn systems. Using numerical simulations, we identify the dominant for- 
mation channels and we predict the expected numbers and characteristics of detectable sys- 
tems, emphasizing how the cluster sources differ from the field population. We explore the 
dependence of formation rates on cluster properties and we explain in particular why the dis- 
tribution of cataclysmic variables has only a weak dependence on cluster density. We also 
discuss the frequency of dwarf nova outbursts in globular clusters and their connection with 
moderately strong white dwarf magnetic fields. We examine the rate of Type la supernovae 
via both single and double degenerate channels in clusters and we argue that those rates may 
contribute to the total SN la rate in elliptical galaxies. Considering coalescing white dwarf bi- 
naries we discuss possible constraints on the common envelope evolution of their progenitors 
and we derive theoretical expectations for gravitational wave detection by LISA. 

Key words: binaries: close - binaries: general - globular clusters: general — stellar dynam- 



1 INTRODUCTION 

From the earliest observations of X-ray binaries in globular clusters 
(GCs) it has been noted that they must be very efficien t sites for the 
production of compact binary systems iClarldri975h . The key to 
the overabundance of compact binaries in clusters, as compared to 
the field, is close stellar encounters. The processes that influence 
the binary population in dense stellar environments include the de- 
struction of wide binaries ("ionization"), hardening of close bina- 
ries, physical collisions, and exchange interactions, through which 
low-mass companions tend to be replaced by more massive partic- 
ipants in the encounter. As a result of these processes, in the dense 
cores of globular clusters, binaries are strongly depleted and their 

i ieriod distribution i s very different from that of a field population 
Ivanova et all2005i) . This effect is stronger for binaries including 
a compact object, like cataclysmic variables (CVs). 

The issue of the dynamical formation of CVs has been exten- 
sively discussed. C onsidering the CV formation via tidal captures, 
ISailvn et alJi 19901) showed that dynamical formation of CVs is not 
expected because more massive donors lead to un s table m ass trans- 
fer. On the other hand, IPi Stefano & Rappaporll il994) predicted 
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the existence of many CVs formed via tidal captures, as many as 
an order of magnitude more than would be predicted by standard 
binary evolution, making CVs a probe of the dynamical processes 
in the cl uster. De tection of CVs in globular clusters proved diffi- 
cult (e.g. IShara et al. 1996) . but a population was detected using 
the Hubble Space Telesc ope iCool et alll995ft . along with a popu- 
lation of "nonflickerers" iCool et aljl998t) which are understood to 
be young helium w hite dwarfs with C/O white dwarf companions 
jHansen et all2003h . 

In the past few years, substantial progress has been made in 
optical identification of Hubble Space Telescope counterparts to 
Chandra X-ray sources in several GCs. Valuable information was 
obtained for populations of CVs, chromospherically active binaries 
and quiescent low-mass X-ray binaries ( qLMXBs) ( Grindlav et^U 



iqu 



'2001VPoo leY et all2002tlEdmonds etltFlaOOSiAlBassa et aj200. 
Heinke et alj|2003h . For the first time we can compare popula 



i 



tions of such binaries in globular clusters (GCs) and in the Galac- 
tic field, and infer their rates of formation and population char- 
acteristics. In particular, 22 CVs have now been identified in 47 
Tuc, allowing identification of several differences between typi- 
cal CVs in globular clusters and CVs in the Galactic field. These 
differences inclu de relatively high X -ray luminosities compared 
to field systems iVerbunt et alJI 19971) : a lack of novae, and of the 
steady, bright blue accretion discs signifying novalike CVs, in GCs 
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iShara & Drissenlll995l) : relatively low frequencies of dwarf nova 
outburst s (DNOs) , the typical identifiers of CVs in the Galactic disc 
^S hara et alll99n) : and a higher ratio of X -ray to optical flux than 
in most field CVs JEdmonds et alj|2003a) . These differences pro- 
duce puzzles: the lack of novae, novalikes, and DNO suggests very 
low mass transfer rates, while the high X-ray luminosities indicate 
moderate mass transfer rates. The X-ray to optical flux ratio sug- 
gests the CVs are DNe, but the lack of DNO argues against this. It 
was suggested that CV discs in GCs are more stable due to a com- 
bination of low mass transfer rates and moderately strong white 
dwarf magnetic moments (Dobrotka et al. 2005). This hints that the 
evolutionary paths of CVs in GCs and in the field are different. 
Comparisons of the numbers of CVs in clusters of different cen- 
tral densities also suppo rts the idea that CV s are produced through 
dynamical interactions JPoolev et alj|2003h . though there is an in- 
dication that CV production may depend more weakly on density 
than t he production of lo w-mass X-ray binaries containing neutron 
stars faeinke et alj2003h . 

This is the first of two papers where we summarize results 
of our studies on compact binary formation in GCs, some pre- 
liminary results of which wer e reported in^yanova et al. (2005); 
llvanova & Rasid i2004l 12003) . In this paper we focus on the for- 
mation of compact binaries with a white dwarf, and in the second 
paper (Paper II) we will describe dynamical formation and evolu- 
tion of binaries with a NS companion. We explore a large spec- 
trum of globular cluster models, where for the first time we take 
into account (i) the mechanism of binary formation through physi- 
cal collisions using results from smoothed-particle hydrodynamics 
(SPH) and (ii) the effect of metallicity on the formation and subse- 
quent evolution of close binaries. In Section 2 we provide a com- 
plete review of the physical processes of formation and destruc- 
tion of mass-transferring WD binaries. In Section 3 we outline the 
methods and assumptions. The major formation channels, and pop- 
ulation characteristics for CVs and AM CVn systems (double WD 
systems where one WD experiences Roche lobe overflow) in differ- 
ent clusters are presented and discussed in Section 4. We conclude 
in the last section by addressing the connection between our results 
and the observations. 



2 MASS-TRANSFERRING WD-BINARIES IN A DENSE 
CLUSTER 

There are several ways to destroy a primordial binary in a globular 
cluster For instance, in a dense region a soft binary will very likely 
be "ionized" (destroyed) as a result of a dynamical encounter. A 
hard binary, in contrast, can be destroyed through a physical colli- 
sion during the encounter. The probability of such an outcome in - 
creases strongly as the binary becomes harder JFregeau et alj2004) . 
In addition to dynamical processes, a primordial binary can be de- 
stroyed through an evolutionary merger or following a SN explo- 
sion. Overall, even if a cluster initially had 100% of its stars in 
binaries initially, the binar y fraction at an age of 10-14 Gyr will 
typically be as low as 10% Cvanova et al 12003) . 

To understand the evolution of a primordial binary in a dense 
environment and the probability of a binary becoming a CV, two 
steps are required: (i) compare the evolutionary time-scales with the 
time-scale of dynamical encounters; (ii) analyze what is the conse- 
quence of an encounter (this depends strongly on the hardness of 
the binary). 

The time-scale for a binary to undergo a strong encounter with 
another single star (the collision time) can be estimated as Tcoii = 



(nEt>oo)~^. Here E is the cross section for an encounter between 
two objects, of masses rrii and mj, with relative velocity at infinity 
Voo and is given as 

S = 7i-d^„^(l + Up/i;^) , (I) 

where dmax is the maximum distance of closest approach that de- 
fines a significant encounter and Vp — 2G{mi + mj)/dmax is 
the velocity at pericenter Assuming that a strong encounter occurs 
when the distance of closest approach is a few times the binary 
separation a, dmax ^ ka with A; ~ 2, we obtain 



r,on = 3.4 x lO" yr fc-^p-^/^M-.'/^n^^^ro' x 



(2) 



1 + 913 



(Mtot + {M)) \ ' 



kPl^Mli^vl, 



Here P^ is the binary period in days, i\ftot is the total binary mass 
in Mq, (M) is the mass of an average single star in A4q, uio — 
?^oo/(10km/s) and 715 = 7i/(10^ pc""^), where n is the stellar 
number density. 

The hardness of a binary system, ry, is defined as 

Gmim2 



V^ 



aa'^{m) 



(3) 



where a is the binary separation, a is the central velocity disper- 
sion, TTii and 7Ti2 are the masses of the binary components, and 
(m) is the average mass of a single star. Binaries that have rj < 1 
are termed soft, and those with 77 > 1 are termed hard. 



2.1 Primordial CVs and AM CVns. 

The typical formation scenario for CVs in the field (low den- 
sity environment) usually involves common envelope (CE) evo- 
lution. In Fig. Q we show parameters of primordial non-eccentric 
binaries that successfully become CVs. To obtain this parameter 
space, we used the binary pop ulation synthesis code StarTrack 
JBelczvnski et all2002ll2005l) ^ We evolved 5 x 10^ binaries con- 
sidering specifically that region of primordial binaries which, ac- 
cording to preliminary lower resolution runs, leads to CV forma- 
tion. Our primary stars have masses between 0.5 Mq and 10 Mq, 
the secondaries have masses according to a flat mass ratio dis- 
tribution with initial periods distributed flatly between I and 10* 
days. For demonstration purposes in Fig.Q we use initially cir- 
cular orbits, because the parameters leading to different formation 
channels can be more clearly distinguished. For our actual cluster 
simulations we use eccentric binaries; in comparison to Fig.Q ec- 
centric primordial binaries can have higher initial periods and still 
produce CVs. Progenitors of CVs with a MS donor are located in 
the left bottom corner, with Alp < 4Mq and logP < 2.5. In 
other cases the donor star is a red giant (RG) or a (subgiant) star 
in the Hertzsprung gap. For primordial binaries located in a small 
but dense area at the left middle part of Fig. \\\ log P ~ 2.7 and 
A/p ~ 1 A'f , a CE does not occur. We note that the lifetime of a 
binary in the CV stage with a RG donor is about 1000 times shorter 
than in the case of a MS donor 

In the core of a GC with core density pc ~ 10* pc~'^, a binary 

^ For the calculatio ns in this study, we hav e used the StarTrack code prior 
to the latest release iBelczvnski et all2005h . However, the most important 
updates and revisions of input physics, in particular the ones important for 
evolution of binaries with white dwarfs, were already incorporated in the 
version we have used. 
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Figure 1. Distribution density of CV progenitors (initial masses of primary 
stars A/p and binary periods P) for non-eccentric binaries in the Galactic 
field, with Z=0.001. The total normalization of CV progenitors is scaled to 
unity, the grey color shows logj^y of the normalized distribution density. The 
thick solid line indicates the binary period where the collision time of the 
binary is equal to the main sequence Hfetime of the primary (using a core 
number density n = 10^ pc~'^, a central velocity dispersion 10 km/s and 
an average object mass of 0.5 Mq). Dash-dotted lines are lines of constant 
binary hardness and dashed lines are lines of constant collision time. 




Figure 2. Distribution density of AM CVn progenitors (initial masses of 
primary stars A/p and binary periods P) for non-eccentiic binaiies in the 
Galactic field, Z=0.001. Notation as for Fig.Q 



with an initial period typical of a CV progenitor will experience a 
dynamical encounter before its primary leaves the MS (see Fig.Q 
where all CV progenitors lie above the line indicating equality of 
Tcou and ta/s). The unaltered primordial channel for CV formation 
is therefore likely to succeed only for binaries that enter the dense 
cluster core after their CE event; the post-CE binary is compact 
enough to avoid an encounter. The contribution of the primordial 
channel depends therefore on the time - before or after the moment 
of CE - when primordial C V binaries will segregate into the central 
dense core. In more detail, an average initial binary in the GC is 
~ O.7M0, which is significantly smaller than the pre-CE mass 
of a primordial CV binary (see Fig. 0. Post-CE primordial CV 
binaries are also heavier than typical binaries in the halo (for which 
the average binary mass is ~ 0.4Af 0). In both cases, primordial 
CV binaries, as heavier objects, will tend to sink toward the cluster 
core on the cluster half-mass relaxation time. 

The situation is similar for the formation of AM CVn systems 
from primordial binaries (see Fig. |2}- In this case, the main for- 
mation channel requires the occurrence of two CE events (see also 
iBelczvnski et alJi2005ah '). and the primordial binary is expected to 
be even wider. However, the second channel, with two stable MT 
stages (at the start of the RG stage of the primary, and when the 
secondary becomes a helium giant), is provided by relatively com- 
pact progenitor binaries. These binaries are expected to evolve in 
the same way in a GC as in the field. 

2.2 Dynamical formation of CVs 

A binary consisting of a MS star and a WD can be formed via sev- 
eral kinds of dynamical encounters: via an exchange interaction, 
via a tidal capture (TC) of a MS by a WD, or via physical collisions 
between a red giant (RG) and a MS star. A fraction of these dynam- 
ically formed MS-WD binary systems will start MT and become a 
CV. In this section we examine in detail the possible channels for 
CV creation. 

The main angular momentum losses in a close MS-WD binary 
occur via magnetic braking (MB) and gravitational wave (GW) 
emission, both of which lead to orbital decay. In eccentric binaries, 
the binary orbital separation will be affected by tides, and the post- 
circularized periastron is larger than the pre-circularized periastron 
(unless tidal synchronization is significant). In Fig.|3|we show the 
maximum initial periods (at the moment of the binary formation) 
of a non-eccentric MS-WD binary that can start MT within 2 Gyr, 
and within 10 Gyr, due only to GW or only to MB (for illustrative 
purposes, we show time-scales for two pr escriptions of magnetic 
braking, one is standard MB according to lRappaport et al] il983h 
(RVJ) and the second is the MB based on dipole-field model ac- 
cording to llvanova & TaamI i2003h (IT03)). A maximum initial pe- 
riod such that a binary is able to start MT without having any other 
encounters is only ~2 days. On Fig.|4]we again show the maxi- 
mum initial periods of binaries that may start MT, but now includ- 
ing all angular momentum losses (GW, MB and tides), and compare 
the cases of non-eccentric and eccentric binaries. On this figure we 
also show the difference in maximum initial period between metal- 
poor and metal-rich GCs. In metal-poor clusters only stars with 
Ad < 0.85 Mq have developed outer convective zones, allowing 
MB and convective tides to operate (lyanova 2006j). This effect can 
potentially be dramatic; for instance, among non-eccentric binaries 
with a MS star of 1 Mq, the range of post-exchange periods that 
leads to C V formation is a factor of 6 larger if the donor has Z=0.02, 
compared to Z=0.001. For eccentric binaries, this ratio is higher, as 
tidal circularization via radiative damping will reduce binary ec- 
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Figure 3. The fate of non-eccentric MS-WD binaries produced by, e.g., 
dynamical encounters, where the primary is a WD of 0.6 Mq. P is the 
post-encounter (or post-CE) orbital period and A/ms is the mass of a MS 
secondary. The short-dashed lines show the binary periods for constant col- 
lision times and the dotted Unes delineate the binaries that will shrink within 
2 and 10 Gyr due to gravitational wave emission. The long-dashed line in- 
dicates the upper period limit for binaries that will begin MT within 2 Gyr 
with the RVJ MB prescription, while the dash-dotted lines indicate those 
that will begin MT within 2 and 10 Gyr with IT03 MB. Below the solid line 
the binary is in contact. 



centricity (and therefore increase the periastron) more effectively 
than GW can shrink the binary orbit. 

A circular binary is most likely to be formed via tidal capture 
(TC), where a post-ca pture circularization is assumed. U sing the 
approach described in IPorteeies Zwart & MeinenI il993h . we can 
estimate the post-capture binary parameters for a MS-WD binary 
(see Fig. 13 where WD mass is assumed to be 0.6 Mq). The upper 
limit here corresponds to the closest approach at which tidal inter- 
actions are strong enough to make a bound system, and the lower 
limit corresponds to the closest approach at which the MS star over- 
fills its Roche lobe by 1/3. We note that the parameter space for 
tidally captured binaries where the MS star does not overfill its 
Roche lobe at the closest approach is very small (see Fig. |3j- We 
note that this is an optimistic estimate, as the captured star can also 
be destroy ed during the chaotic phase of the tidal energy damping 
iMardlinell995.> . Most tidally captured binaries can be brought to 
contact by MB before either the next encounter occurs, or the MS 
star evolves away from the MS. 

An eccentric binary can be formed via an exchange encounter 
or a physical collision; eccentricity can also be increased via the 
cumulative effect of fly-by encounters. For binaries formed through 
MS-RG collisions, the post-exchange binary separation af as well 
as post-exchange eccentricity ef depends on the closest approach p 
^Lombardi et al. 200g) and can be estimated using results of SPH 
simulations. These simulations were done for physical collisions of 
a NS and a RG, and therefore are not straightforwardly applicable 
for the physical collisions of a MS star and RG. We therefore study 
how strongly the choice of the treatment can affect the final results. 
We consider the two following prescriptions: 




Figure 4. The fate of post-encounter MS-WD binaries where the primary 
is a WD of 0.6 A/q, for post-exchange eccentricities and 0.66. P is the 
post-encounter orbital period and A/ms is th^ mass of the MS secondary. 
Thick lines delineate the maximum periods for binaries which will begin 
MT within 2 Gyr. Thin lines of the same type show the period at which that 
binary will begin MT. 



• Using a common-envelope (CE) prescription: 



(Afrg + M,ns)«^ 



■+QCE- 



GM„iMn 



GMrir{Mri 



Af„ 



2ai 



XRi, 



-(4) 



Here A/rg, A/ms and A/wd are the masses of the RG, MS star, and 
RG core that will become a WD, in A/q; Rkg is the RG radius; 
acE is the CE efficiency parameter; and A is the CE parameter 
that connects a star's binding energy with its parameterized form. 
We assume that after a common envelope event the binary is not 
eccentric. 

• Using parameterized results from SPH simulations: 

P 



ef 



Of 



0.1 



SRrg 



3.3(1 



(5) 



(6) 



As the parameterized SPH simulations were done for a limited set 
of mass ratios, we also check the energy balance. When we con- 
sider the case of the second treatment, we choose the minimum bi- 
nary separation from eq. J4} and J6|, as at small masses the extrap- 
olated prescription from SPH simulations can lead to the formation 
of binaries with artificial energy creation. Also, in the case when 
a MS star at the pericenter overfills its Roche lobe, we destroy the 
MS star instead of forming a binary. This is consistent with the re- 
sults of SPH simulations for physical collisions of a RG and a MS 
star (J. Lombardi 2005, priv. communication). 

In Fig.|5]we also show possible binary periods for binaries 
formed via physical collisions with a red giant. Note that it is hard 
to form a relatively close MS-WD binary (one that is able to start 
MT within a few Gyr) with a WD more massive than 0.3 Mq via 
either prescription. Also, in binaries with the mass ratio > 3, Roche 
lobe overflow leads to delayed dynamical instability and a binary 
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Figure 5. Formation of WD-MS binaries via pliysical collisions and tidal 
captures. The hatched area shows binaries formed via TC with a 0.6 Mq 
WD. In the dense hatched area, the MS star did not overflow its Roche lobe 
at the minimum approach during the TC. The dashed lines show binaries 
formed via physical collisions of a MS star and a 0.8 A/q RG, for different 
core masses, using parameterized results of SPH simulations (for illustra- 
tive purposes, we show only the case of the impact parameter to be 0.54 
of the RG radius with corresponding post-collisional eccentricity of 0.7). 
The dotted lines show binaries formed via physical collisions of MS star 
and a 0.8 Mq RG, for different core masses, assuming common envelope 
approach (acB = A = 1). 



merger. This limits the MS star mass to < 0.9 Mq. Therefore, the 
C V progenitors from the channel of physical collisions of RGs and 
MS stars are expected to initially have rather low mass WD ac- 
cretors, and donor star masses < 0.9 Mq. Therefore, metallicity 
variations should not affect this channel strongly. The evolutionary 
stage during He core burning lasts almost the same time as the RG 
branch, however, He core stars of < 2 Mq are a few times more 
compact than at the end of the RG branch and have a larger core 
than during RG evolution. Therefore a collision between a He core 
burning star and a MS star also favors the formation of a WD-MS 
binary that is close enough to become a CV. This channel mainly 
provides binaries with a WD mass at the start of accretion of about 
0.5 Mq (just a bit above the core mass at the time of He core flash). 

On the other hand, there are not many single WDs of such 
small masses present in a GC core. A WD with mass < 0.3 Mq 
cannot (yet) be formed in single star evolution - it must evolve 
via a CE event or a physical collision. A binary containing such a 
WD is very hard and has Tcoii ^ 10 Gyr. If an encounter occurs, 
it is more likely to result in a merger rather than an exchange. We 
therefore expect that most CVs with a low mass WD companion 
will be formed either through a CE event (in a primordial binary or 
in a dynamically formed binary with P ~ 10 — 100 days), or as a 
result of a physical collision, but not via direct exchange encounter. 

A typical binary formed via an exchange encounter has e ~ 
0.7. In order to become a CV within 2 Gyr (or before the next 
encounter), it should have a post-encounter period of a few days 
(see also Fig.|4j. According to energy conservation during an ex- 
change encounter iHegeie et alJll996l) . and assuming that during 
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Figure 6. The fate of post-encounter WD-WD binaries where the primary 
is a WD of 0.6 Mq . P is the post-encounter orbital period and A/wD is the 
mass of the WD secondary. The dashed hnes show the binary periods for 
constant collision times and the dotted lines delineate the binaries that will 
begin MT within 2 Gyr due to GW emission for different post-encounter 
eccentricities. The solid lines show binaries of different eccentricities that 
can be formed through a colUsion between a WD of 0.6 A/q and a RG of 
0.8 Mq ( acEA = 1). 



an exchange encounter the less massive companion is replaced by 
the more massive intruding star, the post-encounter binary separa- 
tion will be larger than pre-encounter. The domain of pre-encounter 
binaries that will be able to form a CV-progenitor binary via only 
exchange encounter is therefore limited to very short period bina- 
ries (with correspondingly long collision times), and these binaries 
are very likely to experience a physical collision rather than an ex- 
change JFregeau et all2004h . 

Let us consider the possibilities for an initially wider dynami- 
cally formed binary than shown on Fig.|4]to evolve toward MT. For 
definiteness, we consider a binary consisting of a 1 Mq MS star and 
a 0.6 A/0 WD with an initial period of 10 days. There are two kinds 
of post-formation dynamical effects that can happen during fly-by 
encounters: (i) binary hardening; (ii) eccentricity pumping. Even if 
each hardening encounter could reduce the orbital separation by as 
much as 50%, the hardening of this binary from 10 days to 1 day (at 
this period MB starts to be efficient) will take about 20 Gyr. In the 
case of eccentricity pumping (assuming no binary energy change), 
the mean time between successive collisions stays at Tcoii ^ 1 Gyr 
and therefore a binary can experience many encounters. If the ac- 
quired eccentricity e ^ 0.95, the binary can shrink through GW 
emission even if its initial period is larger than 10 days. The last 
possibility for such a wide dynamically formed binary to become a 
CV is a CE event that happens in a post-exchange MS-MS binary. 



2.3 Dynamical formation of AM CVns 

For the evolution of WD-WD binaries, we adopt that only GW are 
important as a mechanism of angular momentum loss and neglect 
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the possibility of tidal heating. The maximum possible periods for 
different post-encounter eccentricities are shown on Fig|6| 

Let us first examine a WD- WD binary formation via direct 
exchange. Again, as in the case of MS-WD binaries, a typical ec- 
centricity is e ~ 0.7 and the separation is comparable to the pre- 
exchange separation. The collision time for both pre-encounter and 
post-exchange binaries is so long that both binary hardening and 
exchanges are very rare events. The main difference with MS-WD 
binaries is that post-exchange WD-WD binary periods that will 
allow a binary to evolve to mass transfer (MT) are several times 
smaller for the same eccentricities. Therefore, the exchange chan- 
nel producing a post-exchange binary consisting of two WDs seems 
to be very unlikely. 

A more important channel seems to be the case of an exchange 
encounter that leads to the formation of a MS-WD binary. If the MS 
star is massive enough to become a RG during the cluster lifetime, 
such a binary can evolve through CE and form a close WD-WD 
binary. 

The second important channel is again a physical collision, in- 
volving a single WD with a RG (see Fig.|6| where we show possible 
outcomes of a such a collision). We note that both treatments (pa- 
rameterized SPH results and CE prescription) lead to the formation 
of WD-WD binaries that are roughly equally likely to start the MT. 

We therefore expect that only a post-CE system can become 
an AM CVn, where the post-CE system could be from a primordial 
binary, a post-collision binary, or a dynamically formed binary. 



3 METHODS AND ASSUMPTIONS 

For our numerical simulations of globu lar clusters we use a M onte 
Carlo approach described in detail in llvanova et alj (2003). The 
method couples the binary population synthesis code StarTrack 
^elczynski et al. 2002, 2005), a simple model for the cluster, and 
a small TV-body integrator for accurate trea tment of all relevan t 
dynamical interaction processes (FewBody, iFregeau et all2004) . 
The main update of the code is the treatment of physical collisions 
with a RG, for w hich we now u se the parameterized results of SPH 
simulations from iLombardi et al. 1 2006) as described in §2.2. In 
our code we keep a complete record of all events that happen to 
any cluster star, dynamical (like collisions, tidal captures and ex- 
changes, as well as changes of the binary eccentricity or the binary 
separation after a fly-by encounter), or evolutionary (like common 
envelope events, mass transfers, or SN explosions). This helps to 
analyze the final populations and understand what factors played 
the most significant role in their formation. 

The "standard" cluster model in our simulations has initially 
N = 10^ stars and initial binary fraction of 100%. The distribu- 
tion of initial binary periods is constant in the logarithm between 
contact and 10^ d and the eccentricities are distributed thermally. 
We want to stress here that about 2/3 of these binaries are soft ini- 
tially (the binary fraction provided by only hard binaries gives an 
initial binary fraction of about 20% if the I-D velocity dispersion 
is 10 km/s) and most very tight binaries are destroyed through evo- 
lutionary mergers. Our initial binary fraction is therefore compara- 
ble to the initial binary fractions that are usually used in A'^-body 
codes, where it is assumed for simplicity that very soft binaries will 
not live long as binaries and will only slow down the simulations. 
For more de tailed discussion on t he choice of the primordial binary 
fraction, see llvanova et alJJ2005h . 

For single stars and primarie s we adopted t he broken power 
law initial mass function (IMF) of lKroupal J2002r and a flat mass- 



ratio distribution for secondaries. The initial core mass is 5% of the 
cluster mass and, assuming initial mass segregation, an average ob- 
ject in the core is about twice as massive as an average cluster star. 
At the age of 1 1 Gyr the mass of such a cluster in our simulations 
is ~ 2 X 10^ Mq and is comparable to the mass of typical globular 
clusters in our Galaxy. 

We adopt a core number density ric = 10^ pc~^ (this corre- 
sponds to pc ~ 10**'^ M© pc"^ at the ages of 7-14 Gyr), a half- 
mass relaxation time f,h ~ 1 Gyr and a metallicity Z = 0.001. 
The characteristic velocities are taken as for a King model Wo = 7 
for the cluster of this mass. We take a one-dimensional velocity 
dispersion cti — 10 km/s and an escape velocity from the cluster 
I'osc = 40 km/s. If, after an interaction or SN explosion, an ob- 
ject in the core acquires a velocity higher than the recoil velocity 
^rcc = 30 km/s, an object is moved from the core to the halo. The 
ejection velocity for objects in the halo is i)cj,h = 28 km/s. 

In addition to the "standard" model we also considered cluster 
models with the following modifications: 

• a metal-rich cluster with Z — 0.02 ("metal-rich"); 

• central density ric = lO'* pc~^ ("med-dens") or ric — 
10^ pc"^ ("low-dens") 

• initial binary fraction 50% ("BF05"); 

• RVJ magnetic braking ("fast MB"); 

• treatment of physical collision using a CE prescription ("CE 
coll"); 

• 47 Tuc-type cluster, characterized by a higher density pc — 
lO'^'^ Mq pc"^, higher metallicity Z = 0.0035, ai = 11.5 km/s, 
Wosc = 57 (with the recoil velocity of 52 km/s and t>cj .h =24 km/s) 
and Uh = 3 Gyr ("47 Tuc"). 

"47 Tuc" model describes the GC where currently the largest 
CV population is identified. In order to find a better match with the 
observations, we examined several variations of this model. In par- 
ticular, we considered a model with an initial binary population of 
50% ("47 Tuc-l-BF05"), and a model with an initial binary popula- 
tion is 50%, and the initial core has a smaller mass - 2%, reflecting 
the effect of a longer half-mass relaxation time on the initial popu- 
lation ("47 Tuc-l-SCBF05"). We also examined the sensitivity of the 
final CV production to the CE efficiency parameter, considering the 
case with oceA = 0.1 ("47 Tuc-i-qceA"). 

In order to estimate the effects of dynamics on the popula- 
tion we also ran the same population as in our "standard" model 
(Z=0.001), but without dynamics ("non-dyn"). In order to compare 
to a field population, we considered the population of stars with 
solar metallicity Z=0.02 and with different times of star formation, 
assuming flat star formation rate through last 10 Gyrs ("field"). In 
"non-dyn" model all stars are formed at the zero age, like in GCs. 



4 NUMERICAL RESULTS 
4.1 Formation channels of CVs 

4.1.1 Main formation channels in the "standard" model 

In Fig. we show the formation channels for all CVs that are 
present in a typical cluster (our "standard" model) at the age of 
10 Gyr. Most of these CVs are too dim to be detected, and we con- 
sider separately the population of CVs that can be detectable ac- 
cording to present observational limits, considered specifically for 
the globular cluster 47 Tuc. For the lim iting X-ray luminosity, we 
take I/x > 3 • 10^° ergs s~^ JGrindlav et al. 2001) , and the limiting 
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Figure 7. Main fomiation channels of all CVs that are present in a "standai'd" cluster at the age of 10 Gyr. Notations: numbers in squares represent contributions 
from different channels, small circles contain the labels for each channel. "CE" - common envelope, "Evol. merger" - evolutionary merger, "Evol. ionization" 
- destruction of the binary via SN explosion, "e-f|-" - increase of the eccentricity during fly-by encounters, "Strong binary encounter" - three- or four-body 
encounter with an outcome as exchange, binary destruction ("ionization") or physical collision ("Collision"); "TC" - tidal capture, "Coll RG" - physical 
collision with a RG. 



bolometric luminosity of the donor Ld > O.OGLp,, set by the lim - 
iting magnitude of HST in the cluster core JEdmonds et all2003al) . 

We label channels in the following way: the first character 
indicates the last major dynamical event affecting the binary be- 
fore it becomes a CV in the core, as follows: (1) entering the core 
(primordial binary); (2) companion exchange during a binary en- 
counter; (3) merger during a binary encounter; (4) physical colli- 
sion with a RG during a binary encounter that resulted in a tight 
binary formation with a RG core as a companion; (5) tidal capture; 
(6) physical collision of a single MS star with a single RG. The 
second character indicates a sub-channel by which the binary was 
modified after the last major dynamical event occurred: (a) stands 



for all sub-channels where no strong evolutionary or dynamical 
event occurred; (b) eccentricity of the formed binary was increased 
via binary-binary encounters; (c) common envelope occurred; (d) a 
previous MT episode played the most important role in the orbital 
decay. 

The primordial channel (channel 1) - provides 37% of all 
CVs that are present in the cluster core ( 42% of detectable CVs). 
We call this channel primordial as the binary keeps both its ini- 
tial companions, and no mergers ever occurred to either of them. 
Only 3/4 of CVs formed via this channel are "purely" primordial 
in the sense that they did not experience a significant dynamical 
encounter throughout their life (la, see also Fig.0; most of these 
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"purely" primordial CVs evolved via CE before they entered the 
core. As was predicted in § 2. 1, very few CVs come from the chan- 
nel where CE occurred after a binary entered the core (Ic). 1/5 of 
all primordial CVs would not evolve via CE or start a MT unless 
their eccentricity was increased via fly-by encounters (lb). A small 
fraction of primordial CVs evolved without a CE but with only a 
MT episode on to a MS star (Id), as was described in § 2.1. 

The binary encounters (channel 2, 3 and 4) are responsible 
for the formation of 46% of all CVs, and the same fraction of 
detectable CVs. In most cases the binary that participated in the 
binary encounter was not a primordial binary, but a dynamically 
formed binary. In more than half of cases, a future accretor had 
been a companion in at least 3 different binaries before it acquired 
its final donor. 

The most effective path is the binary exchange channel (chan- 
nel 2) - it provides 32% of all CVs. Within this channel, ~ 40% 
of post-exchange binaries evolved toward the MT without further 
significant dynamical or evolutionary events (2a), in 20% of them 
CE occured (2c) and in 40% of them the MT started as a result of 
the eccentricity pumping during subsequent fly-by encounters (2b). 
This is the most efficient channel for eccentricity pumping. 

Exchange encounters that lead to CV formation typically oc- 
cur between the following participants: (i) a single, relatively heavy 
WD (about 0.7 - 1.4 M©) and a MS-MS binary of total mass 
< 1 Mq ; (ii) a single, relatively massive MS star (about turn-off 
mass) and a MS-WD or WD- WD binary. In the latter case, CE of- 
ten follows the exchange encounter. The number of successful en- 
counters between MS star and WD-WD binary is relatively small, 
and no successful four-body encounter occurred. Nearly all bina- 
ries that proceed via sub-channels 2a or 2b are WD-MS binaries, 
and all binaries in sub-channel 2c are MS-MS binaries after the 
last strong binary encounter. A post-exchange binary typically has 
a heavier WD than a primordial (post-CE) binary has. 

A further 13% of CVs are formed in binaries that experienced 
a physical collision during the last three- or four-body encounter 
- binary coUisional channel (channel 3), while in 1% of cases a 
physical collision with a RG occurred during the encounter and a 
binary with the stripped RG core was formed (channel 4). In the 
evolution of post-collisional binaries the eccentricity change plays 
a smaller role compared to post-exchange binaries; MT is started 
due to the evolutionary angular momentum losses. 

The tidal capture channel (channel 5) contributed very little 
in our standard model. When we looked at all CVs that were formed 
via TC over all GC ages, we find that a typical WD that captured a 
MS star is '^ 1.0 ± 0.2 Mq. In our simulation we allowed a star to 
overfill its Roche lobe radius by up to 1/3 during the tidal capture 
encounter and survive. If all encounters where a MS star overfills 
its Roche lobe lead to the stars' merger, then the contribution of 
tidal captures would be even smaller. 

Finally, the channel of physical collision with RGs (channel 
6) provides 15% of all CVs but much smaller fraction of detectable 
CVs. Eccentricity pumping played a very small role in both TC 
and physical collision channels. Typical participants of a successful 
physical collision (leading to CV formation) are a MS star of 0.3 — 
0.9 Mq and a RG of about 1 — 1.7 A/q with a core around 0.3 Mq 
or a He core burning giant with a core mass around 0.5 Mq. CVs 
formed by this channel are similar to post-CE CVs from primordial 
binaries. We also compared the results of CVs productions in our 
large model with 10^ stars and in the model with three times less 
stars. We noted that, with the increase of the resolution, the total 
number and the number of detectable CVs per unit of the core mass 



is slowly decreasing. Branching ratios between sub-channels within 
a channel can vary slightly, but an overall picture is the same. 
We outline our main findings: 

• Only '^ 25% of CVs were formed in binaries that would be- 
come CVs in the field. 

• In '^ 20% of CVs the main reason for a binary to become 
a CV were fly-by encounters. These CVs cannot be predicted in 
simulations where only strong encounters are taken into account. 

• In ~ 15% of CVs, the WD was formed during dynamical re- 
moval of the RG envelope. As this re moval is not "clean", and about 
0.1 Mq (see iLombardi et al JI2006I) can remain bound to the RG 
stripped core, the characteristics of the WD can differ from those 
formed via a common envelope. 

• 60% of CVs did not evolve via CE, which is the most common 
formation channel for field CVs. 

• Tidal captures did not play a significant role. 



4.1.2 Formation channels in different clusters 

In Table Q we give details on the formation channels for different 
cluster models at the same age of 10 Gyr. We note that these num- 
bers fluctuate with time and are not defined precisely (see more 
below in i|4. 1.3> . however some trends can be identified. We also 
show the numbers of CVs that are formed in a metal-poor environ- 
ment ("non-dyn") and in the field. The definition of "detectable" 
CVs is not very consistent here, as observational limits for field 
CVs are not the same as for GCs (and much more poorly defined), 
but we use the same limits for comparison. It can be seen that dy- 
namics in the "standard" model leads to an increase of the total CV 
production by less than a factor of two. 

In the case of the "metal-rich" model, the turn-off mass at 10 
Gyr is larger than in the "standard" model - there are more massive 
stars in the core; the ratio between the total numbers of CVs in two 
models is roughly the ratio between their turn-off masses at this 
age. 

As was expected, the role of purely primordial CVs (chan- 
nel la) decreases in importance when density increases (see "stan- 
dard", "med-dens" and "low-dens" models), although their abso- 
lute number is about the same for all three models - once a CV 
is formed outside the core, it is hard to destroy it in the core. On 
the other hand, the number of systems that experience CE after en- 
tering the core (Ic) increases as the density decreases, since it is 
easier for pre-CE systems to survive in a less dense environment. 
The production of almost all channels via dynamical encounters de- 
creases with density, except for channel lb, where only non-strong 
encounters are involved. Overall, the total number of CVs in the 
core does not depend strongly on the core density, as the dynamical 
destruction of primordial binaries that would produce CVs, and the 
dynamical production of CVs, compensate each other. 

The "fast MB" model shows the greatest difference with the 
"standard" model in the total number of CVs that are present in the 
cluster core: the "standard" model has about 3 times more CVs, 
both total and detectable, although the number of CVs that are ever 
formed in the core is slightly smaller. The "fast MB" model em- 
ploys the prescription of MB with faster angular momentum loss 
than in the case of the standard model, and therefore the duration 
of the CV stage is shorter. 

The "CE coll" model does not show significant differences 
with our "standard" model. 

The "BF05" model shows that CV formation is reduced 
mainly for primordial CVs and for CVs produced via binary en- 
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Table 1. Formation channels of CVs that are present in the cluster cores of different models at the age of 10 Gyr. 



channel 


la 


lb 


Ic 


2a 


2b 


2c 


3a 


3b 


3c 


4a 


5a 


6a 


Total 


Detec 


standard 


0.271 


0.077 


0.013 


0.135 


0.129 


0.052 


0.090 


0.006 


0.039 


0.013 


0.026 


0.148 


209 


47 


metal-rich 


0.204 


0.056 


0.031 


0.148 


0.143 


0.046 


0.051 


0.015 


0.015 


0.020 


0.031 


0.230 


265 


16 


med-dens 


0.327 


0.253 


0.167 


0.111 


0.012 


0.056 


0.031 


0.019 


0.006 


0.006 


0.000 


0.006 


193 


35 


low-dens 


0.404 


0.066 


0.456 


0.037 


0.000 


0.007 


0.015 


0.000 


0.015 


0.000 


0.000 


0.000 


156 


26 


fast MB 


0.190 


0.103 


0.017 


0.086 


0.190 


0.172 


0.034 


0.017 


0.000 


0.017 


0.052 


0.103 


79 


15 


CE coll 


0.212 


0.106 


0.006 


0.159 


0.194 


0.041 


0.041 


0.029 


0.029 


0.006 


0.000 


0.176 


230 


47 


BF05 


0.206 


0.119 


0.024 


0.135 


0.135 


0.056 


0.040 


0.032 


0.024 


0.024 


0.008 


0.175 


162 


36 


47Tuc 


0.135 


0.094 


0.000 


0.250 


0.146 


0.073 


0.042 


0.042 


0.021 


0.031 


0.000 


0.156 


275 


37 


47 Tuc+BFOS 


0.143 


0.057 


0.014 


0.171 


0.143 


0.029 


0.014 


0.043 


0.029 


0.000 


0.043 


0.300 


190 


35 


47 TUC-I-SCBF05 


0.071 


0.114 


0.000 


0.100 


0.114 


0.057 


0.014 


0.014 


0.014 


0.014 


0.029 


0.443 


237 


27 


47 Tuc-l-acE A 


0.170 


0.057 


0.011 


0.182 


0.125 


0.045 


0.034 


0.011 


0.011 


0.023 


0.023 


0.307 


253 


40 


non-dyn 


























124 


16 


field 


























117 


3 



Notations for channels - see text in S I4. 1.1 l and also Fig.Q "Total" is the number of CVs and "Detec" is the number of detectable CVs, both numbers are 

scaled per 50 000 A/q stellar population mass in the core. 




5xl0» 10" 

Age [yr] 
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Figure 8. Formation of CVs via different channels. The solid line shows the 
case of "standard" model, dotted line shows "metal-rich". 



Figure 9. Appearance of CVs (time that MT starts) formed via different 
channels. Notations as in Fig.lsl 



counters. The number of CVs produced via physical collisions be- 
tween single stars is about the same. 

The results for "47 Tuc" are due to a mixture of several condi- 
tions: the higher core number density favors dynamical formation, 
and the higher metallicity gives a wider mass range over which MB 
operates on the donor star. Variation of initial conditions, such as a 
smaller binary fraction, does not lead to significant differences ex- 
cept that the relative role of binaries becomes smaller than the role 
of physical collisions. Some decrease in the number of detectable 
CVs occurs when we start with a smaller initial core. 



4.1.3 Formation channels at different cluster ages 

In Fig.|8|we show the formation rate of CVs via different channels 
throughout the life of a cluster (indicating the time when a dynam- 



ical event occurred, or, for primordial binaries, the time when a CE 
happened). In Fig.|9|we show the rate of appearance of CVs (start 
of mass transfer) as a function of time. 

The primordial channel produces most of its CVs at the be- 
ginning of the cluster evolution, though the appearance of primor- 
dial CVs is distributed flatly in time. This contrasts with binary en- 
counter channels, where the formation occurs rather flatly in time, 
but the rate of CV appearance grows after 7 Gyr. A similar delay in 
the appearance can be seen for CVs formed via tidal captures. 

In Fig. llOl we show the number of CVs that stop MT for dif- 
ferent reasons: (1) end of MT due to the star's contraction (usually 
occurs with CVs where a donor is a RG) or evolutionary merger of 
the binary; (2) explosion of the WD as la SN or sub-Chandrasekhar 
SN explosion, or accretion induced collapse (AIC); (3) end of MT 
due to a strong dynamical encounter. Most CVs stop MT due to an 



© 0000 RAS, MNRAS 000, 000-000 



10 A'^ Ivanova et al. 



Table 2. Formation channels of CVs that are present in the "standard" cluster core at different ages. 



channel 


la 


lb 


Ic 


2a 


2b 


2c 


3a 


3b 


3c 


4a 


5a 


6a 


Total 


Det 


IGyr 


0.250 


0.000 


0.000 


0.000 


0.250 


0.250 


0.000 


0.000 


0.000 


0.000 


0.000 


0.250 


15 


7 


2Gyr 


0.111 


0.074 


0.037 


0.259 


0.111 


0.185 


0.000 


0.000 


0.037 


0.000 


0.000 


0.185 


84 


37 


3Gyr 


0.185 


0.046 


0.015 


0.154 


0.169 


0.123 


0.015 


0.000 


0.015 


0.031 


0.015 


0.231 


170 


62 


4Gyr 


0.200 


0.078 


0.022 


0.167 


0.133 


0.133 


0.022 


0.000 


0.011 


0.022 


0.011 


0.200 


203 


76 


5Gyr 


0.190 


0.076 


0.029 


0.162 


0.086 


0.095 


0.057 


0.010 


0.019 


0.019 


0.019 


0.229 


210 


62 


6Gyr 


0.214 


0.077 


0.026 


0.162 


0.103 


0.060 


0.051 


0.009 


0.026 


0.026 


0.017 


0.231 


211 


56 


7Gyr 


0.268 


0.049 


0.016 


0.163 


0.114 


0.041 


0.049 


0.008 


0.024 


0.024 


0.016 


0.220 


204 


43 


8Gyr 


0.284 


0.061 


0.014 


0.128 


0.122 


0.068 


0.054 


0.007 


0.041 


0.020 


0.014 


0.182 


228 


44 


9Gyr 


0.268 


0.067 


0.013 


0.128 


0.134 


0.060 


0.074 


0.007 


0.047 


0.020 


0.013 


0.161 


214 


40 


lOGyr 


0.271 


0.077 


0.013 


0.135 


0.129 


0.052 


0.090 


0.006 


0.039 


0.013 


0.026 


0.148 


209 


47 


llGyr 


0.253 


0.084 


0.006 


0.157 


0.139 


0.054 


0.078 


0.006 


0.042 


0.012 


0.030 


0.139 


212 


51 


12Gyr 


0.209 


0.088 


0.005 


0.165 


0.181 


0.049 


0.066 


0.016 


0.033 


0.005 


0.027 


0.154 


221 


41 


13Gyr 


0.203 


0.091 


0.005 


0.188 


0.162 


0.056 


0.066 


0.030 


0.030 


0.010 


0.025 


0.132 


228 


34 


14Gyr 


0.199 


0.131 


0.000 


0.184 


0.155 


0.039 


0.068 


0.044 


0.029 


0.010 


0.019 


0.121 


229 


43 



Notations are as in TableHI 
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Figure 10. Destruction of CVs. Notations as in Fig.|8| 



Figure 11. The mass-distribution of hydrogen accreting WDs. The hatched 
area coiTesponds to dynamically formed binaries; the solid filled area to 
systems formed directly from primordial binaries. The top panel shows the 
case with no dynamics (different ages), the second panel from the top shows 
the compiled field case, the third panel shows the core of the cluster in the 
"standard" model and the bottom panel shows the cluster core in the model 
"47 Tuc". 



evolutionary reason, while the number of SN explosions is also rel- 
atively high and is comparable to the number of CVs destroyed by 
dynamical e ncounters (for more det ail on how SN la are calculated, 
see §4.4 and lBelczvnski et alfcOO^ . 

In Table|2|we show a detailed representation of CV formation 
by different channels in the "standard" model at different cluster 
ages. We note a peak in the number of detectable CVs at the age 
of 4 Gyr, and that the total number of CVs increases steadily until 
the age of 8 Gyrs and then stays constant. The weight of different 
channels in the relative numbers of appearing CVs does not change 
dramatically during the cluster evolution. 



4.2 Population characteristics of CVs 

4.2.1 Periods and masses 

On Fig. llll and ll2l we show mass and period distributions for "stan- 
dard" and "47 Tuc" cluster models, as well as for the population 
evolved without dynamics. A remnant of the primordial popula- 
tion of CVs in cluster cores follows the distribution of primordial 
CVs evolved without dynamics at the same age. However the dis- 
tribution of dynamically formed CVs shows signs of younger (non- 
dynamical) CV populations, and also is more populated at the high 
mass end. For 47 Tuc we also show the periods of the identified 
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Figure 12. The period-distribution of hydrogen accreting WDs. The 
hatched area corresponds to dynamically formed binaries; the solid filled 
area to systems formed directly from primordial binaries. The top panel 
shows the case with no dynamics (different ages), the second panel from 
the top shows the compiled field case, the third panel shows the core of the 
cluster in the "standard" model and the bottom panel shows the cluster core 
in the model "47 Tuc". Solid triangles indicates the periods of CVs that are 
identified in 47 Tuc from observations. 



Figure 13. Simulated distribution of orbital periods versus 0.5-2.5 keV X- 
ray luminosity for "47 Tuc" model (at the age of 1 1 Gyr). Stars are CVs 
that could be detected, while open circles are CVs that cannot be detected 
(generally due to the optical faintness of the donor). The size of the symbol 
corresponds to the WD mass (the largest is for WDs more massive than 1 
Mq, the smallest symbol is for WDs less massive than 0.6 A/q, and the 
medium symbol is for WDs with masses between 0.6 and 1 A/0 ). 



CVs, which have a distribution consistent with the period distribu- 
tion of the "detectable" CVs in our simulations. 



4.2.2 MT rates and X-ray luminosities 

We compare the distribution of orbital periods versus 0.5-2.5 keV 
X-ray luminosity for our simulation of 47 Tuc and observations of 
CVs in globular clusters (see Fig. 1 131 and Fig. I14t ^. Observation- 
ally, few CVs have measured orbital periods. For real CVs with un- 
known periods, only the X-ray luminosity is shown. To obtain the 
X-ray luminosity of simulated CVs for comparison to observations, 
we use the accretion model from P atterson & Raymond ( 1985) for 
0.5-4 keV and scale the luminosity to 0.5-2.5 keV assuming a flat 
energy distribution within the band: 



Lx(0.5 - 2.5 keV) = 0.066 



2i?wd 



(7) 



where M is the mass transfer rate and _Rwd is the radius of the WD. 
The simulations show reasonable agreement with the observations. 

However, this picture so far does not explain the rare occur- 
rence of DNOs in globular clusters CVs, in comparison to the field. 
Therefore, other properties of these systems must be explored to 
identify the distinguishing characteristics. 

In Fig. llSl we show MT rates in CVs in "standard" model, in 
the "47 Tuc" model and for the field. Although MT rates in our GC 
simulations do not exceed 10~^ Mo/yr, and field CVs can have 



Note that our "47 Tuc" represents a cluster similar to 47 Tuc but 5 times 



less massive 
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higher MT rates, this result may be due to small number statis- 
tics. We therefore do not find significant differences between MT 
rates in GC CVs and field CVs. For our GC CVs, all MT rates are 
such that the accretion disc is partially ionized and unstable, in ac- 
cordance with the disc instability model. We adopt the viscosity 
parameters ahot = 0.1 and Ocoid = 0.01, for hot and stable disc 
state s respectively. Therefo re, all our CVs should produce DNOs. 

iDobrotka et alJi2005l) proposed that a combination of low MT 
rates and moderately strong magnetic fields can explain the absence 
of DNOs in GCs. Lower mass transfer rates would lead to a rarer 
occurrence of DNOs, and strong magnetic fields would lead to trun- 
cation of the inner disc keeping the disc in the cold stable state. As 
we do not find systematically lower MT rates for our GC CVs, (in 
fact our MT rates are t wo orders of magnitude higher than found in 
iDobrotka et alJ i2005h ). we estimated the minimum magnetic field 
required to suppress D NOs (see Fig. ll6> . using the criterion from 
jDobrotkaetalJJ2005h ): 



Bsupp > 5.7xlO''G 



M 



10- 






O.Oli?, 



rCw 



.(8) 



We find that Bsupp is a slowly increasing function of the WD mass 
and, for most of the WDs with masses below 1 Mq, is even below 
lO'' G. WDs with B < lO'' G are not regarded as highly magnetic. 
It can also be seen that a 10^ G field is enough to prevent DNOs 
in all CVs with the WDs less massive than 1 . 1 AIq and the field of 
10* G is strong enough to stop DNOs for WDs of all masses. 



4.3 Formation Channels of AM CVns 

In Table 3 we show the main formation channels for AM CVns 
that occur in different clusters, classifying channels in an analo- 
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Table 3. Formation channels of AM CVns that are present in the "standard" cluster core at 10 Gyr. 



channel 


la 


lb 


Ic 


2a 


2b 


2c 


3a 


3b 


3c 


4a 


5a 


6a 


Total 


standard 


0.385 


0.068 


0.000 


0.051 


0.090 


0.030 


0.013 


0.013 


0.064 


0.064 


0.009 


0.265 


316 


metal-rich 


0.379 


0.089 


0.000 


0.031 


0.071 


0.022 


0.018 


0.004 


0.009 


0.009 


0.000 


0.348 


303 


med-dens 


0.732 


0.141 


0.056 


0.010 


0.010 


0.015 


0.005 


0.000 


0.000 


0.000 


0.000 


0.076 


236 


low-dens 


0.680 


0.085 


0.178 


0.008 


0.004 


0.036 


0.000 


0.000 


0.000 


0.000 


0.000 


0.178 


283 


fast MB 


0.450 


0.077 


0.000 


0.036 


0.068 


0.045 


0.009 


0.005 


0.027 


0.027 


0.000 


0.266 


303 


CEcoU 


0.474 


0.123 


0.000 


0.035 


0.082 


0.023 


0.006 


0.012 


0.018 


0.018 


0.012 


0.211 


231 


BF05 


0.271 


0.090 


0.005 


0.053 


0.080 


0.016 


0.011 


0.032 


0.011 


0.011 


0.000 


0.404 


242 


47Tuc 


0.246 


0.070 


0.000 


0.026 


0.070 


0.035 


0.035 


0.018 


0.009 


0.009 


0.000 


0.491 


327 


47 Tuc-l-BFOS 


0.174 


0.110 


0.000 


0.018 


0.101 


0.037 


0.000 


0.009 


0.009 


0.009 


0.009 


0.514 


296 


47 TUC-I-SCBF05 


0.217 


0.101 


0.000 


0.029 


0.043 


0.043 


0.000 


0.000 


0.014 


0.014 


0.000 


0.551 


233 


47 Tuc-HacE A 


0.252 


0.078 


0.000 


0.000 


0.117 


0.029 


0.019 


0.010 


0.058 


0.058 


0.000 


0.427 


296 


non-dyn 


























169 


field 


























110 



Notations for channels is the same as for CVs (see text in S I4. 1.1 l and also Fig.0. "Total" is the number of AM CVns, the number is scaled per 50 000 Mq 

stellar population mass in the core. 



gous manner to Sec. 4.1.1. The main differences with the forma- 
tion of CVs are: (i) post-CE channels (including primordial) are 
more important; (ii) the role of physical collisions for AM CVns 
is increased; (iii) binary encounters played a much less significant 
role. 

In most physical collisions the participants are a RG of 0.9 — 
2.2 Mq (the core mass is 0.2 - 0.3 Mq) and a WD > 0.6 A/q. In 
< 20% of physical collisions the participants are a He core burning 
giant of 1.7 — 2.2 Mq (the collision in this case leads to the forma- 
tion of a He star of 0.4 — 0.55 Mq) and a heavier WD, generally 
> 0.9-1.1Mq. 

Overall, this channel provides AM CVns with accretors of 
masses from 0.8 Af© to 1.4 Af© at a cluster age of 10 Gyr. 
The field population of AM CVns would have accretors with 
masses typically between 0.65 Mq and 1.0 M q (Fig. 1 171 see also 
iNelemans et all200ltlBelczvnski et all2005d) . The peak in the ac- 
cretor mass distribution is shifted from ~ 0.7 Mq in the field pop- 
ulation to 0.9 Mq in the cluster population. 



4.4 Explosive events 

As the mass distribution of accreting WDs is shifted towards higher 
masses compared to the field population, it is important to check 
what effect this has on the rates of: (i) Type la supernova (SN la, 
here we mean the single degenerate channel only); (ii) double WD 
mergers (those for which the total mass ^ A'/ch — 1.4 A/q); (iii) 
sub-Chandrasekhar supemovae; and (iv) accretion induced collapse 
(AIC). 

The type of event that occurs depends on the mass and com- 
position of the white dwarf and the rate of mass transfer. If this 
is a carbon-oxygen (CO) WD and experiences stable accretion, 
it will accumulate mass until it reaches the Chandrasekhar limit 
and then explodes as a type la supernova. In the case of specific 
MT r ates the accretion lead s to the accumulation of He in the 
shell teato & Hachisiill999l) . If sufficient mass is accumulated, it 
will lead to the ignition of the CO or ONeMg core and disrup- 
tion ofjheWD as a sub-Chandrasekhar-mass Type la supernova 
(seelTaanJl98(llLivne & GlasneJl99il:IWooslev & WeaveJl994L 
iGarcia-Senz et alll999tllvanova & Taaml2004^ . In the case of ac- 
cretion on to ONeMg WDs, upon reaching the Chandrasekhar limit 
the WD will undergo accretion-induced collapse and form a NS. 



Table 4. Rate of explosive events. 



event 


SNIa 


supraCh 


subCh 


AIC 


NSdd 


NSaic 


standard 


2.0 


9.3 


7.3 


3.21 


118 


79.2 


metal-rich 


3.48 


8.2 


11.7 


2.73 


117 


91.4 


med-dens 


0.48 


2.4 


8.0 


1.45 


89 


64.7 


low-dens 


0.00 


2.9 


9.8 


0.72 


80 


61.3 


fast MB 


3.48 


8.2 


11.7 


2.73 


117 


91.4 


CE coll 


0.74 


7.7 


8.9 


2.97 


117 


85.1 


BF05 


1.56 


3.6 


4.2 


2.45 


70 


45.5 


47Tuc 


1.22 


4.9 


9.3 


2.44 


88 


70.3 


47 TUC-HBF05 


1.54 


2.2 


6.4 


0.66 


65 


42.1 


47 TUC-HSCBF05 


0.66 


2.4 


6.2 


0.88 


50 


36.3 


47 Tuc-HQcE-^ 


2.20 


3.2 


8.8 


1.71 


87 


59.6 


non-dyn 


0.00 


0.7 


3.2 


0.00 


70 


69.4 


field 


0.50 


7.7 


15.1 


1.58 


122.17 


22.8 



"SN la" is the number of Type la SN (single-degenerate channel only), 
"supraCh" is the number of double WD mergers where the total mass is 
more than 1.4 AfQ, "subCh" is the number of sub-Chandrasekhar nuclear 
runaways, "NS" is the number of NSs that can be potentially formed via 
double-degenerate (NSqij) or AIC (NSaic) channels until the age of 10 
Gyrs. For cluster models, rates and numbers are given per Gyr per 200,000 
A/0 total cluster mass and are averaged for the ages of 8-12 Gyr. For non- 
dynamical models, rates and numbers are given for the age of 10 Gyr and 
for the field model rates and numbers are given after 10 Gyr of continuous 
star formation. 



If the donor is another WD and the mass transfer is not sta- 
ble, the mass of the merger product can exceed the Chandrasekhar 
limit - these so-called supra-Chandrasekhar mergers could lead ei- 
ther to a Type la supernova (double-degenerate channel), or to a 
merger induced collapse of the remnant to form a N S and perhaps 
a millisecond radio pulsar iChen & Leonar j|l993h . It was argued 
as well that in the latter case and, if one of the WD is magnetic, 
such mergers will lead to magnetar formation. Such objects may be 
responsible for the production of the giant flares emitted by soft 7- 
repeaters, which can be identified with early type galaxies. These 
flares may contribute to a fraction of the observed s hort duration 
burst population at higher redshift JLevan et al.l2006h . 

If NSs are bom with natal kicks, most of them will be ejected 
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Figure 14. Observed distribution of orbital periods vs. 0.5-2.5 keV X- 
ray luminosity for CVs in globular clusters. Known CVs in globular clus- 
ters, with known X-ray luminosities but without known orbital periods, 
are plotted at P=0 hours or less. CVs that have undergone recorded dwarf 
nova outbursts are indicated with circles, CVs with strong He II A4686 
emission (suggesting strong magnetic fields, see text) with stars (some 
are both). Three objects with uncertain periods or CV status (W34 in 47 
Tuc might be a millisecond pulsar) are marked as open triangles; other 



CVs are indicated with filled triangles. Data incl ude 23 CVs from 47 
Tuc ( Edmonds et al."2003a'b'). 8 from NGC 6397 "Edmonds et al.""l99^ 
iGrindla ' 



iGrindlav et al. 2001: Kaluznv & Thompson 2003; Shara et al. 20O 



from NGC 6752 IPooley et al JI200 2 ; Bailvn et al. 1996), two from M22 
iPietrukow icz et alj| 2005h. two from M15 iHannikainen et alJ|2005t). one 
from M5 iHakala et alll997HNeill et all2002l) . one from M4 ifiassa et alj 
l2004l) . and one from M55 JKaluznv et all20 



from the shallow cluster potential, leavin g few NSs to exp lain the 
observed number of millisecond pulsars JPfahl et alj|2002f) . In the 
case of accretion or merger induced co llapse, the NSs are likely t o 
be formed without a significant kick JPodsiadlowski et al]l2004^ ■ 
and this can relieve the NS "retention problem". If double WD 
mergers do not lead to collapse, they must contribute to the rate of 
Type la sup emovae, with pote ntial cosmological implications (for 
a review see[Leibundgut^200 11) . 

The production of supra-C handr asekhar mergers in our galaxy 
was discussed in lHurlev et alJ i2002h and was estimated to be 2.6 
per year in the Galactic disc (this is 8.6 per cluster per Gyr in our 
units). Several free parameters can have strong effects on this re- 
sult, such as the common envelope prescription, the initial mass 
function, or the adopted star formation his tory. There are also dif- 
ferences between our models and those of iHurlev et alj i2002ft : (i) 
their cut-off mass for WD binaries is at the initial mass of 0.8 Mq ; 
(ii) they choose acE = 3; (iii) they adopted continuous star forma- 
tion through 15 Gyr (c.f. ours 10 Gyr); and (iv) our model for accre- 
tionon to WDs is more up-to-date (for details, see Belczynski et al. 
l2005h . Overall, we find that our form ation rates for the fie ld are not 
significantly different from those of iHurlev et alj yOO^) (see Ta- 
ble |4}. We also find that our rates are smaller if the star formation 
is taken not as flat, but with one (or several) star formation bursts 
that ended several Gyrs ago. 



Figure IS. MT rates in CVs. Hatched areas indicate detectable CVs. For our 
no-dynamics model the detectable CV histogram is increased by a factor of 
10, while for "47 Tuc" it is increased by a factor of 3. "47 Tuc" is shown at 
the age of 1 1 Gyr 



The enhanced production rate of double WD mergers in dense 
stellar clusters was first discussed in detail by IShara & HurlevI 
J2002r . who applied this to open clusters. They found that the 
supra-Chandrasekhar WD merger rate can be increased by an or- 
der of magnitude (although their statistics were based on only a 
few events). We did not find such an increase compared to the field 
population, where star formation is continuing, though we found 
some increase compared to the case without dynamical interactions 
(see Table |4j. However, we note that our total number of supra- 
Chandrasekhar WD mergers is large. In fact, if indeed all those 
mergers lead to formation of NSs, and those NSs are retained by 
the cluster, then this channel provides about 6% of all NSs ever 
created. The NSs thus created become comparable in numbers to 
the NSs that were born with natal kicks and retained. The produc- 
tion of NSs via this channel can be reduced by reducing the effi- 
ciency of the common envelope. In this case, more binaries will 
merge during the CE phase and less supra-Chandrasekhar mergers 
will occur. We found however that even the reduction of qceA to 
0.1 led only to a moderate decrease of the "current" (at about 10 
Gyr) production rate of supra-Chandrasekhar mergers, while their 
total production is only a bit smaller (see different models for 47 
Tuc in the Table|4j). In addition, the production of NSs via AIC is 
comparable to the production of NSs via merger induced collapse, 
and therefore also appears to be a significant source of NSs in GCs. 
The question of how many NSs can be created via different chan- 
nels in GCs is very important, and will be addressed in more detail 
in Paper II. 

We estimate the contribution of SN la produced in GCs to to- 
tal galactic SN rates. Assuming that ~ 3 x 10^ Mq is contained 
in galactic GCs (~ 150 galactic GCs), we find that the single- 
degenerate channel from GCs can provide at best 1 SN per ~ 10^ yr 
per galaxy, and the contribution of GCs is only several times higher 
if the double-degenerate channel (supra-Chandrasekhar) also re- 
sults in SN. In spiral galaxies the rate of SN la is 0.04-0.1 per cen- 
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Figure 16. The minimum magnetic field Bsupp required to prevent dwarf 
nova outbursts for CVs in our simulations. Open circles show CVs in the 
"standard" model, filled triangles show CVs in the "metal-rich" model, and 
stars indicate CVs in the "47 Tuc" model. "47 Tuc" is shown at the age of 
1 1 Gyr, other models are shown at the age of 10 Gyr. 
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Figure 17. The mass-distribution of helium accreting WDs. The hatched 
area corresponds to dynamically formed binaries; the solid filled area to 
systems formed directly from primordial binaries. The top panel shows the 
case with no dynamics (different ages), the second panel from the top shows 
the compiled field case, the third panel shows the core of the cluster in the 
"standard" model and the bottom panel shows the cluster core in the model 
"47 Tuc". 
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Figure 18. Distribution density (averaged over time) of LISA binaries, in 
the space of binary periods and chirp masses, from our "standard" model 
(integrated over cluster ages from 9 to 13 Gyr). 



tury per galaxy jMannucci et all2005l) . and therefore the contribu- 
tion of GCs is not important. However, GCs can play a larger role in 
elliptical galaxies, where no star formation is going on and the rate 
of SN la provided by the field population is smaller. In addition, in 
ellipticals the specific frequency of GCs per galaxy luminosity unit 
is significantly higher than in spirals (up to 8 .6 compared to 0.5 
in Milky Way, see e.g. iKim & Fabbiandl2004) . Also, it has been 
shown that the observational SN la rate consists of two compo- 
nents - a prompt piece that is proportional to the star formation 
rate, a nd an extended piece that is pr oportional to the total stellar 
mass JScannapieco & Bildstenl 120051) . This is consistent with the 
behavior of the formation rates of both single degenerate and dou- 
ble degenerate channels in GCs, which peak during the first Gyr of 
the cluster evolution and have a flat distribution at later times. We 
therefore propose that GCs can increase the theoretically predicted 
rates of SN la in elliptical galaxies. 



4.5 Coalescing compact binaries as LISA sources 



AM CVns discussed in i|4.3l as well as most double white dwarf 
mergers discussed in !|4.4l (with the exclusion of a small fraction of 
mergers that are results of physical collision during hard binary en- 
counters) are coalescing binaries driven by gravitational radiation. 
Prior to their mergers, or before and during the MT, they can be de- 
tectable as gravitational wave sources by LISA. Their detectability 
is significantly enh anced when their orbital periods become sma ller 
than about 2000 s JBenacouista et alfeOoHJStroeer et aj2005h . so 
their signals can be distinguished from the background noise pro- 
duced by Galactic binaries. As the positional accuracy of LISA will 
be much greater for binaries with such short periods, these sources 
can be associated with specific globular clusters in our Galaxy. 

From our simulations we find that at any given moment, a typ- 
ical cluster of 200,000 Mq will contain 10 LISA sources on aver- 
age, and at least 3 LISA binaries at any given moment; during 1 
Gyr a cluster forms 180 LISA systems. A massive cluster like 47 
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Tuc, with mass r^ lu iviq, 
any given moment, and 40 LISA binaries on average. At an age of 
10 Gyrs such a cluster can produce as many as 3-15 NS-WD LISA 
binaries per Gyr (a typical cluster produces 1-3 NS-WD coalescing 
binaries per Gyr). 

With the total mass in GCs of about 3 x 10^ Mq, as many as 
1500 LISA binaries can be present in all Galactic GCs (this number 
will decrease if the CE efficiency is smaller). The most optimistic 
upper limit for the galactic formation rate of NS-WD binaries in 
GCs is several hundred per Gyr. The lifetime of a NS-WD binary 
in the LISA-band during the MT is ~ 10* yr. The time prior to the 
onset o f MT depends on the binary eccentricity, but is usually much 
shorter Jivanova et alJ2005h . With our predicted formation rates, as 
many as 10-50 NS-WD LISA binaries can be detected in GCs. This 
number is probably too optimistic, as fewer than 10 ultra-compact 
X-ray binaries (mass-transferring NS-WD binaries with orbital pe- 
riods less than an hour) have been identified in GCs, implying that 
the formation rate of LISA-sources should be smaller. (One expla- 
nation could be that many globular clusters are less dense than our 
"standard" model and NS-WD formation is thus less frequent.) We 
shall address the issue of the formation rates of binaries with NSs 
in more detail in Paper II. Here we will only note that the LISA bi- 
naries will spend spend most of their time in the LISA band among 
their MT tracks, with chirp masses < O.25M0 (see Fig. llSt . 



5 DISCUSSION 

With our simulations we predict that the formation rates of C Vs and 
AM CVns in globular clusters are not very different from those in 
the field population. The numbers of CVs and AM CVns per mass 
unit are comparable to numbers in the field if the whole cluster 
population is considered, and they are only 2-3 times larger in the 
core than in the field. Dynamical formation is responsible only for 
60%-70% of CVs in the core. This fraction decreases as the density 
decreases, and the role of primordial CVs becomes more impor- 
tant. We rule out tidal captures as an effective mechanism for CV 
formation in GCs unless the rate of TCs is significantly underesti- 
mated. Instead we propose that the population of GC CVs reflects 
a combination of primordial CVs, CVs in post-exchange binaries, 
and products of physical collisions of MS stars with RGs. There are 
also primordial CVs which are located in the halo and have never 
entered the core. The GC core density variation indeed does not 
play a very large role, in contrast to the case o f NS binaries, where 
almost all systems are formed dynamically ilvanova et alJl2005l) 
and whose numbers have a str ong dependence on the cluster col- 
lision rate JPoolev et alJl2003B . We expect to have one detectable 
CV per 1000 Mq in the core of a typical cluster and about one de- 
tectable CV per 1000 - 2000 Mq in a 47 Tuc type cluster. Thus 
we predict 35-40 CVs in the core of 47 Tuc, in quite reasonable 
agreement with observations, where 22 CVs in 47 Tuc have been 
identified ( Edmonds et al. 2003aj). Even better agreement between 
our simulations and the observed number of CVs can be obtained 
if we assume that the initial core mass in 47 Tuc is smaller than 5% 
of the cluster. 

Although the formation rates do not differ strongly, we found 
significant differences in the populations, and note that these differ- 
ences may have observational consequences. Indeed, cataclysmic 
variables in globular clusters have an unusual array of character- 
istics, which make them difficult to classify as members of the 
standard classes of CVs recognized in the galaxy. Their X-ray 
luminosities seem to be rather high, compared with CVs in the 



field iVerbunt et alll997n . They exhibit dwarf nova outbursts only 
rarely, compared to well-studied dwarf novae: only 1 dwarf nova 
was found in 47 Tuc by Shara et al. ( 1996) in a survey which would 
have detected 1/3 of known dwarf no vae if they were located in 
47 Tuc, while lEdmondsetalJ J2003J) identify 22 firm CVs in 47 
Tuc. Finally, the X-ray to optical flux ratios of CVs in globular 
clusters are relatively high, comparable to those of dwarf novae 
JEdmonds et al. 2003bi) . 

One solution to this problem was the suggestion that CVs in 
globular clusters tend to be primaril y mag netic in nature, compared 
to CVs in the field Icrindlav et alJll995h . Magnetic CVs have no 
discs (AM Her or polar CVs), or truncated discs (DQ Her CVs or 
intermediate polars, IPs), because of the effect of the WD mag- 
netic field. As a result, the disc instability is nonexistent or sup- 
pressed. Magnetic CVs are believed to produce X-rays through an 
accretion shock above the polar cap, producing high X-ray lumi- 
nosities, while nonmagnetic CVs produce a n optically thick bound- 
ary la yer, saturating their X-ray emission JPatterson & Ravmond 
Il985l) . Strong He II A 4 686 lines were observ ed in the spectra of 
three CVs in NGC 6397 JEdmonds et alll999^ ■ indicating a strong 
source of FUV radiation. This FUV radiation could indicate either 
evidence for an intermediate polar interpretation, or a very high 
mass transfer rate; the second interpretation is favored for the FU V - 
bright, 27-hour period CV AK09 in 47 Tuc JKnigge et al1l2003h . 
Another argument in favor of the intermediate polar interpretation 
is the excess Nh (in addition to that expect ed along the line o f 
sight) observed towards many CVs in 47 Tuc jHeinke et al 120051) . 
Excess Nh in CVs that are not observed at high inclinations has 
been considered a signature of the accretion curtains observed in 
the magnetic systems known as intermediate polars. 

However, only two globular cluster CVs have shown clear 
evidence of magne tic fields in their X-ray lightcurves so fa r (X9 
and XIO in 47 Tuc lGrindlav et alJl200lLlHeinke et alj|2005h . This 
may not mean that these systems are not magnetic, since the num- 
ber of X-ray photons detected from globular clusters is generally 
small (compared with nearby, well-studied CVs). In addition, it 
has been suggested that the accretion in SW Sex and VY Scl CVs 
is governed by the WD magnetic field, without evidence of pul- 
sations ( Rodriguez-Gil et al. 2001 ; Hameurv & Lasota 200.3 ). An- 
other problem for the magnetic interpretation is that IPs tend to be 
optically brighter than typical CVs in globular clusters, which have 
lower ra tios of X-ray to optica l flux more typical of dwarf novae 
than IPs JEdmonds et all2003d) . A final problem is the observation 
of dwarf nova outbursts in tw o of the three CV s in NGC 6397 pos- 
sessing strong He II emission JShara et all20 05:). A proposed reso- 
lution to these problems is a combination of a low mass transfer rate 
(which will reduce the optical brightness and incre ase the X-ray to 
optical flux ratio) with an interm ediate polar nature JEdmonds et alj 
l2003d) . iDobrotka et alj J2005h calculated the dwarf nova recur- 
rence times for CV discs with various mass transfer rates and WD 
magnetic moments, and found a parameter space that fulfilled the 
requirements of globular cluster CVs. Left unanswered was why 
globular cluster CVs might tend to have stronger magnetic fields 
than field systems. 

Our work provides a possible answer to this question. Glob- 
ular cluster dynamics has a strong effect on the composition 
of the binaries that form CVs, tending to place more massive 
WDs into binaries that will become CVs (Fig. 11). Increasing 
the mass of WDs in CVs increases the energy that can be ex- 
tracted at a given mass transfer rate, thus increasing the X-ray 
luminosity and X-ray to optical flux ratio of the CVs. This ef- 
fect is complementary to the effects of higher magnetic fields. 
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However, higher mass WDs a lso have a higher probabihty of 
showing strong magnetic fields (Vennes 1 999l : lLiebert et all2003l : 
[Wickramasinghe & Ferrario 2005 ). Thus the dynamical origin of 
WDs in globular cluster CVs may be responsible for the obser- 
vational peculiarities of globular cluster CVs; their relatively high 
X-ray luminosities and X-ray to optical flux ratios, and their low 
rates of dwarf nova outbursts. 

The tendency for higher mass white dwarf accretors in GC 
CVs in comparison to the field also affects the production of the 
superhump phenomenon. This behavior results from the precession 
of the outer disc due to the excitation of resonances within the disc 
caused by the 3: 1 commensurabil ity of motions in the disc w ith the 
companion's orbital period (see IWhitehurst & King 1 19911) . Such 
systems are characterized by mass ratios (of donor to accretor) of 
less than 0.25-0.33. Systems of this type in the field are rarely ob- 
served at orbital periods above the period gap, but the higher white 
dwarf masses of CVs in GCs would increase their likelihood in 
GCs. 

We examined also several other consequences of having a dy- 
namically modified population of close binaries including WDs. 
In particular, considering supra-Chandrasekhar mergers, we found 
that too many NSs may be formed if these mergers lead to merger- 
induced collapse. We suggest that either this mechanism does not 
lead to NS formation, or the CE efficiency is overestimated. By our 
estimates, GCs do not contribute strongly to the SN la rates in spi- 
ral galaxies, however they may significantly increase these rates in 
elliptical galaxies. We have also shown that GCs can be excellent 
sites for LISA observations since many GCs will contain several 
LISA sources at any given moment, although most of those sys- 
tems will have low chirp masses. 
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